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ABSTRAGT 

Recent suggestions of a “photon underproduction crisis” (Kollmeier et al. 2014) have generated con¬ 
cern over the intensity and spectrum of ionizing photons in the metagalactic ultraviolet background 
(UVB). The balance of hydrogen photoionization and recombination determines the opacity of the 
low-redshift intergalactic medium (IGM). We calibrate the hydrogen photoionization rate (Ph) by 
comparing Hubble Space Telescope spectroscopic surveys of the low-redshift column density distribu¬ 
tion of H I absorbers and the observed (z < 0.4) mean Lya flux decrement, Da — (0.014)(1 J- z)^'^, 
to new cosmological simulations. The distribution, /(iVnuz) = (Pj\f/d{\ogNui)dz, is consistent with 
an increased UVB that includes contributions from both quasars and galaxies. Our recommended fit, 
rH(z) = (4.6 X 10“^^ s“^)(l -I- z)^-^ for 0 < z < 0.47, corresponds to unidirectional LyG photon flux 
$0 « 5700 cm“^ s“^ at z = 0. This flux agrees with observed IGM metal ionization ratios (G III/G IV 
and Si Ill/Si IV) and suggests a 25-30% contribution of Lya absorbers to the cosmic baryon inven¬ 
tory. The primary uncertainties in the low-redshift UVB are the contribution from massive stars in 
galaxies and the LyG escape fraction (/esc), a highly directional quantity that is difficult to constrain 
statistically. We suggest that both quasars and low-mass starburst galaxies are important contribu¬ 
tors to the ionizing UVB at z < 2. Their additional ionizing flux would resolve any crisis in photon 
underproduction. 

Subject headings: cosmological parameters — ultraviolet: galaxies — observations — intergalactic 
medium — quasars: absorption lines 


1. INTRODUCTION 

One the most important but poorest known parame¬ 
ters in studies of the intergalactic medium (IGM) and cir- 
cumgalactic medium (CGM) is the intensity of the meta¬ 
galactic ionizing ultraviolet background (UVB). This un¬ 
certainty is not surprising, considering the strong ab¬ 
sorption of extreme ultraviolet (EUV) photons at wave¬ 
lengths A < 911.753 A by the interstellar medium (ISM). 
Lacking direct measurements of the metagalactic radia¬ 
tion field in the Lyman continuum (LyG), astronomers 
rely on indirect probes of photoionizing radiation and 
theoretical estimate^ based on the presumed ionizing 
sources and cosmological radiative transfer. An accu¬ 
rate characterization of this UVB is necessary for model¬ 
ing the thermal and ionization conditions (Bolton et al. 
2014; Boera et al. 2014) of the IGM and CGM, for a cen¬ 
sus of IGM baryon content (Shull et al. 2012b), and for 
the ionization corrections needed to derive metallicities 
(Shull et al. 2014; Werk et al. 2014). Por photoionization 
of hydrogen, helium, and many spectroscopically acces¬ 
sible ions of heavy elements (C, N, O, Ne, Mg, Si, S, Fe) 
the most important spectral range extends from the hy- 
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^ UVB models include calculations by Haardt Sz Madau (1996, 
2001, 2012), denoted HM96, HMOl and HM12, Shull et al. (1999), 
and Faucher-Giguere et al. (2009). HM05 refers to the higher-flux 
spectrum in a 2005 August update to the Haardt Sc Madau (2001) 
“Quasars-|-Galaxies” spectra, provided for inclusion in the pho¬ 
toionization code Cloudy. The hydrogen photoionization rates, 
Fh, ill these models vary by up to a factor of six (see Table 1). 


drogen Lyman edge (energies E > 13.6 eV) through the 
He II continuum {E > 54.4 eV) and into the extreme 
ultraviolet (EUV) and soft X-ray (100-1000 eV). 

The claim of a “photon underproduction crisis” 
(Kollmeier et al. 2014, hereafter denoted K14) generated 
considerable anxiety in the extragalactic astronomy com¬ 
munity and prompted a re-examination of the assump¬ 
tions in modeling the EUV background at low redshift. 
These authors compared predictions from smoothed- 
particle hydrodynamic (SPH) simulations to observed 
properties of low-redshift H I (Lya) absorbers. The bal¬ 
ance of hydrogen photoionization and radiative recombi¬ 
nation rates sets the opacity of the low-redshift IGM, 
which can be measured from the distribution of Lya 
absorbers in redshift (z) and H I column density, A/ji 
(cm“^). This allows one to infer the unidirectional ion¬ 
izing photon flux $0 (cm“^ s“^) and specific intensity U 
(erg cm“^ s“^ Hz“^ sr“^) of LyC radiation. K14 found 
that the hydrogen photoionization rate (Th) required to 
match their simulated distribution of H I absorbers was 
five times larger than the value predicted in a recent cal¬ 
culation of the UVB by HM12. They found better agree¬ 
ment with previous estimates of a higher UVB (Shull 
et al. 1999; HMOl) that included significant contribu¬ 
tions from galaxies. A large part of the UVB discrepancy 
can be traced to the anomalously low values at z < 2 of 
the LyG escape fraction, /esc = (1-8 x 10“'^)(1 -I- x)^'^, 
adopted by HM12 for galaxies. 

In this paper, we undertake a careful examination of 
this UVB discrepancy, comparing observations of low-z 


2 


TABLE 1 

Hydrogen Photoionization Rates^ 


Model Reference 

Tnfz = 0) 

ThIz = 0.25) 

Haardt Sz Madau (1996) 

4.14 x 

8.63 x 10-13 

Shull et al. (1999) 

6.3 X lO-i"* 

15 x 10-13 

Haardt & Madau (2001) 

10.3 x 

20.6 x 10-13 

Haardt &; Madau (2005) 

13.5 x 

29.1 x 10-13 

Haardt & Madau (2012) 

2.28 x 

5.89 x 10-13 

Faucher-Giguere et al. (2009) 

3.84 x 

7.28 x 10-13 

Shull et al. (2015) 

4.6 x 10-13 

12 x 10-13 


^ Hydrogen photoionization rates Eh (s computed at redshifts 
2 = 0 and 2 = 0.25 by various theoretical models. 


UVB extrapolated to low redshift was under-estimated 
by HM12 by a factor of approximately 2-3. 

2. CONSTRAINING THE UV BACKGROUND 

2 . 1 . Definitions and Measurements of the Ionizing 
Radiation Field 

For an isotropic radiation field of specific intensity 
Ii,, the normally incident photon flux per frequency is 
(irl^/hv) into an angle-averaged, forward-directed effec¬ 
tive solid angle of tt steradians. The isotropic photon 
flux striking an atom or ion is Air^y/hv), and the hy¬ 
drogen photoionization rate follows by integrating this 
photon flux times the photoionization cross section over 
frequency from threshold {vf) to oo. 


Lya absorber to new cosmological simulation^ Differ¬ 
ences among previous theoretical estimates for the UVB 
may arise from modeling of the evolution of sources of 
EUV radiation from galaxies and quasars. We use H I 
data from our recent surveys of low-redshift Lya ab¬ 
sorbers with UV spectrographs on the Hubble Space Tele¬ 
scope {HST): the STIS (Space Telescope Imaging Spec¬ 
trograph) survey of 746 Lya absorbers (Tilton et al. 
2012) and the COS (Cosmic Origins Spectrograph) sur¬ 
vey of over 2600 Lya absorbers (Danforth et al. 2015, 
updated in July 2015). These surveys give consis¬ 
tent results for the bivariate distribution, /(Vhi,z) = 
dfN/d{\ogN}ii)dz, of absorbers in redshift {z < 0.47) 
and H I column density (12.5 < loglVni < 15.5). 

In Section 2, we describe current estimates of the 
UVB and low-redshift hydrogen photoionization rates 
and compare our HST surveys of intergalactic H I column 
densities with new grid-code simulations of the IGM. We 
calibrate the UVB through its influence on the distri¬ 
bution of H I column densities and the flux decrement 
{Da) from Lya line blanketing. In Section 3, we justify 
boosting Th by a factor of 2-3 above HM12 values, to¬ 
gether with other parameters that characterize the UVB: 
the specific intensity /q (erg cm“^ s“^ Hz“^ sr“^) at 
the Lyman limit and the integrated LyC photon flux 
$0 (cm“^ s“^). For 0 < z < 0.47, our revised ion¬ 
ization rate is Fh ~ (4.6 x s“^)(l -I- z)'^'^. At 

z = 0, this corresponds to $0 « 5700 cm“^ s“^ and 
Iq = 1.7 X 10“^^ erg cm“^ s“^ Hz“^ sr“^. These param¬ 
eters are in agreement with photoionization modeling in 
our recent IGM survey (Shull et al. 2014) of metal-ion 
ratios (C III/C IV and Si Ill/Si IV). This increase in 
the UVB is comparable to previous calculations (Shull 
et al. 1999; HMOl, HM05) and likely arises from impre¬ 
cise modeling of sources of EUV radiation (galaxies and 
quasars). The larger UVB would be consistent with an 
increased contribution of galaxies with LyC escape frac¬ 
tions, /esc ~ 0.05 at z < 2, in contrast to the very low 
values /esc < 10“^ adopted by HM12. A recent analy¬ 
sis (Khaire & Srianand 2015) of the effects of a revised 
QSO luminosity function (Groom et al. 2009; Palanque- 
Delabrouille et al. 2013) increased the low-redshift UVB 
by a factor of two. Both of these results suggest that the 

^ Our simulations are made with the N-body hydrodynamic grid 
code Enzo (http://enzo-project.org). For an overview of the code 
see Bryan et al. (2(J14j. Uur previous applications of Enzo to IGM 
astrophysics appear in Smith et al. (2011), Shull et al. (2012a), and 
Shull et al. (2012b). 
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Here, we approximate the frequency dependence of spe¬ 
cific intensity and photoionization cross section by power 
laws, Ii, = Iq{v/vq)~°‘ and Uy = ao{iy, where 
(Jo = 6.30 X 10“^® cm^ and a « 1.4 for AGN (Shull 
et al. 2012c; Stevans et al. 2014). The integrated unidi¬ 
rectional flux of ionizing photons is then 


4*0 = 



dv 


tt/o 

ha 


( 2 ) 


which is related to the density of hydrogen-ionizing pho¬ 
tons by $0 = n-f{c/A) for an isotropic radiation field. We 
then have the relations among parameters: 


Lh = 4(7o4>o 
Anioao 


= (8.06 X 10-1^ s-^)$4 (3) 


rH = 


h{a + 3) 


a J- 3 

= (2.71 X 10-14 s-i)/_23 


( 4 ) 


where we normalize the incident flux of ioniz¬ 
ing photons and specific intensity to characteris¬ 
tic values, $0 = (10^ cm-^s-i )$4 and Iq = 

(10“^® erg cm”^ s”! Hz”! sr-i )/_23 at the hydrogen Ly¬ 
man limit {hvo = 13.60 eV). In an HST/COS survey of 
IGM metallicity at z < 0.4, Shull et al. (2014) found that 
$4 « 1 and /_23 « 3 gave reasonable fits to the observed 
ratios of adjacent ionization states of carbon and silicon, 
(Si Ill/Si IV) = 0.67;[( f^, (G III/G IV) = 0.70;((;2o, and 

their sum, (Hem + Hciv)/(f 2 siiii + Hsiiv) = 4.9tei. 

Over the past 20 years, numerous papers have es¬ 
timated the ionizing background and photoionization 
rate. Table 1 lists Fh for several models, with val¬ 
ues at z = 0 ranging from (2.28 — 13.5) x lO-^^ s”!. 
The HM12 rate, Fh = 2.28 x lO-^^ s”!, corresponds 
to a one-sided ionizing flux <I>o = [rH(a -|- 3)/4 ctoQ:] ~ 
2630 cm”^ s”! and specific intensity I -23 ~ 0.823 for 
the radio-quiet AGN spectral index (a = 1.57) assumed 
by HM12. These fluxes are lower by a factor of three 
compared to the z = 0 metagalactic radiation fields 
from AGN and galaxies calculated by Shull et al. (1999), 
/agn = X 10-23 and Jcai = l.ltj;® x 10-23, re¬ 

spectively. Adding these two values with propagated 
errors gives a total intensity and hydrogen ionization 
rate of Rot = 2.4tJ;g x 10-23 erg cm -2 s-^ Hz-^ sr-^ 
and Fh = 6 . 0 / 2 ’^ x IQ-^^ s-^. The difference between 
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TABLE 2 

Column Density Distribution^(STIS Survey) 


(log Ahi) 

Range in log A^hi 

abs 



13.0 

(12.9-13.1) 

115 

4.637 

124l« 

13.2 

(13.1-13.3) 

115 

5.003 


13.4 

(13.3-13.5) 

106 

5.216 

102lj« 

13.6 

(13.5-13.7) 

75 

5.327 

70l| 

13.8 

(13.7-13.9) 

73 

5.341 

68l® 

14.0 

(13.9-14.1) 

50 

5.360 

47l® 

14.2 

(14.1-14.3) 

35 

5.379 

3316 

14.4 

(14.3-14.5) 

26 

5.382 

2416 

14.6 

(14.5-14.7) 

18 

5.382 

17l^ 


^ For 613 low-redshift Lyo absorbers taken from HST/STIS 
survey (Tilton et al. 2012), the columns show: (1) mean column 
density (A^hi in cm“^) with bin width AlogA^Hi = 0.2; (2) 
bin range in logA^ni; (3) number of Lyo absorbers (A/^bs) in 
bin; (4) total redshift pathlength A^eff at each Aijji; (4) bivari¬ 
ate distribution of absorbers in column density and redshift, 
/(iVniiz) = d?Af/d{\ogNiii)<iz, computed as (A/’abs/0.2 A^gff) 
and plotted in Figures 2-7. 

these radiation fields appears to be the contribution from 
galaxies. The HM12 background adopts a negligible 
UVB from galaxies, whereas the background models of 
Shull et al. (1999), HMOl, and HM05 have comparable 
intensities from galaxies, owing to higher assumed LyC 
escape fractions. 

2.2. HST Observations of the HI Column Density 
Distribution 

Previous ultraviolet spectroscopic surveys of low- 
redshift Lya absorbers estimated their contribution to 
the baryon census through the distribution of H I col¬ 
umn densities in the diffuse Lya forest (e.g., Penton 
et al. 2000, 2004; Lehner et al. 2007; Danforth & Shull 
2008). Our recent HST surveys of the low-redshift IGM 
were more extensive, obtaining 746 Lya absorbers with 
STIS (Tilton et al. 2012) and 2577 Lya absorbers with 
COS (Danforth et al. 2015). The COS survey used the 
medium-resolution far-UV gratings (Green et al. 2012) 
with coverage between 1135-1460 A (G130M) and 1390- 
1795 A (G160M); a few spectra extended slightly out¬ 
side these boundaries. Our COS survey probed 82 AGN 
sight lines with cumulative pathlength Az = 21.7 with 
H I column densities A^hi (in cm“^) between 12.5 < 
logiVni < 17. For this paper, we analyze a “uniform 
redshift-limited sample” of Lya absorbers at z < 0.47. 
As discussed in these survey papers, the Lya statistics 
are consistent with 24-30% of the baryons residing in the 
Lya forest and partial Lyman-limit systems. 

Table 2 summarizes the data selected from the STIS 
survey, with 613 Lya absorbers at z < 0.4 over the range 
12.9 < log TVhi < 14.7, reprocessed as described in Dan¬ 
forth et al. (2015). Table 3 gives similar results from 
the COS survey at z < 0.47, with 2074 absorbers be¬ 
tween 12.6 < logiVni < 15.2. For each bin in logA^ni, 
we list the number of absorbers {Af^hs) and effective red¬ 
shift pathlength (Azes) over which our Lya survey is 
sensitive. From these data, we derive the bivariate dis¬ 
tribution, d'^Af/d{logN}ii)dz, of absorbers in H 1 col¬ 
umn density and redshift by dividing Afahs by AzeS and 


TABLE 3 

Column Density Distribution^ (COS Survey) 


(log Ahi) 

Range in logA^Hi 

abs 



12.7 

(12.6-12.8) 

194 

3.67 

2601)66° 

12.9 

(12.8-13.0) 

292 

9.79 

isoltg 

13.1 

(13.0-13.2) 

383 

16.46 

120 ± 10 

13.3 

(13.2-13.4) 

368 

18.78 

98 ± 5 

13.5 

(13.4-13.6) 

267 

19.23 

69l)^ 

13.7 

(13.6-13.8) 

146 

12.03 

61 ± 5 

13.9 

(13.8-14.0) 

123 

12.25 

50 ± 5 

14.1 

(14.0-14.2) 

96 

12.96 

37 ± 4 

14.3 

(14.2-14.4) 

75 

13.75 

271" 

14.5 

(14.4-14.6) 

47 

14.17 

17tl 

14.7 

(14.6-14.8) 

45 

14.24 

wfl 

14.9 

(14.8-15.0) 

19 

14.25 

fi 7+1-9 
^+-1.5 

15.1 

(15.0-15.2) 

19 

14.25 

n 7+1-9 
^+-1.5 


^ Statistics of low-redshift Lyo absorbers (12.6 < log Ahi < 15.2) 
taken from the HST/COS survey (Danforth et al. 2015, updated 
July 2015) using the “uniform sub-sample” of 2074 H I absorbers. 
The columns show: (1) mean column density (A^hi in cm“^) with 
bin width AlogAni = 0.2; (2) bin range in logApii; (3) number 
of Lyo absorbers (AAbs) in bin; (4) total redshift pathlength Azeff 
at each Ahi; (4) bivariate distribution of absorbers in column 
density and redshift, /(Ahi,2) = /d{log Niii)dz, computed as 
(A/'abs/0.2 Azeff) and plotted in Figures 2-7. 


AlogAni = 0.2. An accurate calculation of this dis¬ 
tribution function requires both good absorber counting 
statistics and knowledge of the pathlength AzeS covered 
in the survey for a given column density. As in Dan¬ 
forth & Shull (2008), we compute Azeff corresponding to 
the 4(7 minimum Lya equivalent width as a function of 
wavelength in each spectrum. Asymmetric error bars 
arise from statistical uncertainties in N^hs and Az^g, 
computed for each bin using the formalism of Gehrels 
(1986). At low column densities, the errors are domi¬ 
nated by Azeff, while at high column densities the errors 
are dominated by the small values of fVabs in bins at 
loglVni > 14.6. 

An important effect of the IGM is the flux decrement. 
Da , produced by Lya forest line blanketing of the contin¬ 
uum flux. Previous measurements of Da were reported 
by Kirkman et al. (2007) using low-resolution UV obser¬ 
vations with the HST Faint Object Spectrograph (FOS) 
toward 74 AGN. Their measured decrements were fitted 
to a power law, Da{z) = (0.016)(l-|-z)^ °^, over the range 
0 < z < 1.6 with considerable scatter. Based on our 
moderate-resolution HST/GOS survey. Figure 1 shows 
the fraction of light removed by Lya absorption from the 
continuum for 0 < z < 0.47. These statistics were found 
by summing the observed-frame equivalent widths of all 
identified (> 4 (t) Lya absorbers in the Danforth et al. 
(2015) catalog within a given redshift range, normalized 
by the effective redshift pathlength probed by the survey 
for absorbers of that strength. The uncertainty, (7 _d^, is 
the quadratic sum of measured equivalent-width uncer¬ 
tainties normalized in the same manner. The errors are 
small at lower redshifts (z < 0.35) where most of the sur¬ 
veyed Lya absorbers are found. However, they remain 
small compared to Da even at higher redshifts where 
Azeff is smaller. Figure 1 shows our results for bins of 
width Az = 0.01 and Az = 0.05 to illustrate the vari- 
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ance on small scales. The scatter in Da is significantly 
larger than and cosmic variance is significant, par¬ 
ticularly at z > 0.35. These fluctuations are large in the 
Az = 0.01 bins, which are comparable to the mean ab¬ 
sorber separations for the observed Lya absorption-line 
frequencies dN/dz « 100 — 200 for logiVni = 12.7—13.3. 

Figure 1 shows a small dip in Da between 0.2 < 
0 < 0.3. This redshift band corresponds to Lya ab¬ 
sorbers at 1459 A to 1580 A, observed primarily by the 
COS/G160M grating (1390-1795 A) since G130M ends 
at 1460 A. The dip appears in several ranges of column 
density, 13 < logAni <14 and 14 < logAni <15. The 
first range is expected to dominate line blanketing in the 
Lya forest for line profiles with Doppler velocity param¬ 
eters 6 « 25 — 30 km s“^. The dip also remains when 
we offset the redshift bins. Thus, we believe the dip to 
be real, although we do not have a plausible physical 
explanation for its presence. 

We fitted the mean values of the HST observations 
with Az = 0.05 bins to the power-law form Da{z) = 
(0.014 ± 0.001)(1 -b z)2-2±0-2. For z < 0.2, these ob¬ 
servations are in reasonable agreement with the lower- 
resolution HST/FOS results of Kirkman et al. (2007), 
although we find steeper redshift evolution. By post¬ 
processing the H I absorption-line profiles in our simula¬ 
tions (Sections 2.3 and 2A) over the range 0 < z < 0.2, 
we find mean decrements at (z) = 0.1 of Da = 0.0116, 
0.0097, and 0.0332 for the HMOl, HM05, and HM12 ra¬ 
diation fields, respectively. At (z) = 0.1, these UVBs 
correspond to hydrogen photoionization rates Fh = 
3.54 X 10-14 s-i (HM12), 1.38 x lO-^^ s"! (HMOl), 
and 1.86 x lO”!^ s”! (HM05). Our simulations find that 
Da oc Fgi’ i’. In their SPH simulations, K14 found decre¬ 
ments at z « 0.1 of Da — 0.024 (HMOl) and Da = 0.050 
(HM12), higher than our simulations by factors of 1.9 and 
2.9, respectively. Similar offsets are seen in the distribu¬ 
tion of H I column densities. 

Our observations at z < 0.47 with STIS (Tilton et al. 
2012) and GOS (Danforth et al. 2015) yield consistent 
values for /(TVhi, z), the bivariate distribution of Lya ab¬ 
sorbers in column density and redshift. A least-squares 
fit of the GOS data over the range 12.7 < logiVni < 
15.2 had the form, /(Ahi,z) = d?M/d{\ogNYii)dz « 
( 167 )iVj) 3 ° ®^4=0 02 ^ where we define the dimensionless col¬ 
umn density Ni^ = [Ahi/ 101^ cm”^]. We can use this 
line frequency to make an analytic estimate of the Lya 
line blanketing and flux decrement, 

= (^)/(iV,^)rf(logiV) 

2 \ pN2 

-) N^r'^dNi^) 

In the last expression of Equation (5), we assume un¬ 
saturated Lya lines, with equivalent widths W\/\ = 
{-Ke^/mec){Nf\/c), where W\ and A = 1215.67 A are 
defined in the rest-frame and / = 0.4164 is the Lya 
oscillator strength. We adopt limits logiVi = 12.5 
and logA ’2 = 14.0 for the integral, which we denote 
I{Ni,N 2 ) « 4.49. The estimated decrement Da = 
(0.00323)/(Ai,iV 2 ) « 0.0145, close to our observations 
at z « 0. Because line frequency /(Ahi,z) oc (1 -b z)^'^ 


/7re^/AA /167 x 10^^ cm 
V mec2 ) 2.303 



Fig. 1. — Flux decrement from the Danforth et al. (2015) 
low-z IGM survey with COS (updated in July 2015) showing 
the fraction of flux removed from the continuum by Lya 
absorbers for redshifts 0 < z < 0.47. The decrement Da{z) 
is calculated from the normalized, summed,observed-frame 
equivalent widths of Lya absorbers in a given redshift bin, di¬ 
vided by the effective pathlength (Azeu) probed by the survey 
at each redshift. Red open circles show Da in bins of width 
Az = 0.01. Filled black circles show Da in bins of Az = 0.05 
and fitted to Da{z) = (0.014)(1 -b 2)^'^4=0-2 line). 

Blue dotted line shows the fit, Da{z) ~ (0.016)(1 -b z)^'°^, 
found by Kirkman et al. (2007) in lower-resolution HST/FOS 
measurements. The small drop in Da at z « 0.2 — 0.3 
appears over several column-density ranges and remains 
when we shift redshift bins. 


out to z « 0.5 (Danforth et al. 2015), we expect Da to 
have the same dependence, as seen in Figure 1. This cal¬ 
culation shows that for the observed steep distribution of 
H I column densities, most of the line blanketing occurs 
from moderate-strength Lya absorbers at Nm « N 2 , be¬ 
ginning to saturate and appear on the flat portion of the 
curve of growth. Owing to line saturation, the stronger 
but rarer Lya absorbers add a small amount to this es¬ 
timate. 

2.3. Cosmological Simulations of HI in the low-z IGM 

Our simulations of the H I absorber distributions were 
run using the Eulerian N-body -b hydrodynamics code 
Enzo (Bryan et al. 2014). The N-body dynamics of the 
dark matter particles in our simulations were calculated 
using a particle-mesh solver (Hockney & Eastwood 1988), 
and the hydrodynamic equations were solved using a 
direct-Eulerian piecewise parabolic method (Golella & 
Woodward 1984; Bryan et al. 1995) with a Harten-Lax- 
van Leer-Gontact Riemann solver (Toro et al. 1994). Our 
modifications of this code and its applications to the IGM 
are discussed in Smith et al. (2011), where we conducted 
tests of convergence, examined various feedback schemes, 
and compared the results to IGM thermal phases, O VI 
absorbers, and star-formation histories using box sizes 
of 25/1-4 Mpc and 50h-^ Mpc and grids of 256A 384A 
512^, 768^, and 1024^ cells. We note that the Smith et al. 
(2011) simulations used the HM96 background, rather 
than HMOl as stated in that papei[^ This code has a 

^ Because the HM96 background is 2.5 times lower than HMOl, 
the LycK forest was stronger in the Smith et al. (2011) simulations. 
However, as we discuss in Section 2.4, the HM96 background is 
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substantial IGM heritage, with applications to O VI ab¬ 
sorbers (Smith et al. 2011), the baryon census at z < 0.4 
(Shull et al. 2012b), clumping factors and critical star- 
formation rates during reionization (Shull et al. 2012a), 
and synthetic absorption-line spectra for comparison of 
simulations with HST observations (Egan et al. 2014). 

Our current simulations were initialized at a redshift 
z = 99 and run to z = 0, using the WMAP-9 maxi¬ 
mum likelihood concordance values (Hinshaw et al. 2013) 
with = 0.282, Oa = 0.718, = 0.046, Ho = 

69.7 km s“^ Mpc“^, erg = 0.817, and rig = 0.965 to 
create the initial conditions. The dark matter density 
power spectrum used the Eisenstein & Hu (1999) trans¬ 
fer function. Three independent realizations of the ini¬ 
tial conditions were created in order to constrain the ef¬ 
fects of cosmic variance due to the finite box size. Two 
of the realizations were created using Music (Hahn & 
Abel 2011) with second-order Lagrangian perturbation 
theory. The third realization was created using the ini¬ 
tial condition generator packaged with Enzo. The basic 
simulations discussed in this paper were run on static, 
uniform grids with a box size of 50h~^ comoving Mpe 
and 512^ and 768^ cells. We also analyzed a previous 
simulation with 1536^ cells, produced by Britton Smith 
on the XSEDE supercomputer and discussed in our study 
of IGM clumping factors (Shull et al. 2012a), as well the 
(1024^, 50/i“^ Mpe) simulations of Smith et al. (2011). 
Because those 1024^ models were run with an older UVB 
from Haardt & Madau (1996), we chose not to compare 
them to our 512^ and 768^ runs, which used more recent 
backgrounds (HMOl, HM05, HM12). 

In view of the different predictions of our code for 
the distribution of low-redshift Lya absorbers, it is ap¬ 
propriate to compare the resolutiorj^ of our grid-code 
simulations to the SPH simulations run by K14. With 
a co-moving box size of 50/i“^ Mpe, our models have 
spatial resolutions (cell sizes) of 97.6/i“^, 65.1/i“^, and 
32.6h~^ kpc for grids of 512^, 768^, and 1536^, respec¬ 
tively. The corresponding dark-matter mass resolutions 
are rridm = 6.1 x 10^/i"! Mq (5123), g ^ 10^/i"! Mq 
(7683), 2.2 X 10® Mq (15363). For comparison, 

the simulations analyzed by K14 used the SPH code of 
Dave et al. (2010), also in a 50h~^ Mpe box with 576® 
dark-matter (and baryon) particles. The characteristic 
spatial resolution in their SPH method is the interparticle 
distance, 50/i“^ Mpc/576 « 86 .8h~^ kpc. Under gravi¬ 
tational evolution, the SPH particles are concentrated in 
regions with baryon overdensity A;, = pb/pb > 1, and 

the interparticle distance is reduced by a factor of Aj,' . 
Although K14 do not quote a mass resolution, their simu¬ 
lation (50 /i“ 3 Mpe, 5763) ^ould have rudm ~ 5 x lO^/i”^ 
Mq, scaling from values quoted in the (5123, xoO/i”^ 
Mpe) SPH calculation of Tepper-Garcia et al. (2012). 

close to the value we recommend in our current study. Conse¬ 
quently, our 2011 simulations are in reasonable agreement with 
the HST observations of the low-redshift Lya forest. 

The original study by Dave et al. (2010) used simulations in a 
Mpe box, and most of their studies were run with 384^ ps.r- 
ticles. To investigate effects of numerical resolution and box size, 
they also employed a simulation with 512^ particles in a 96/i“^ Mpe 
box. Contemporaneous SPH studies (Tepper-Garcia et al. 2012) of 
the IGM were conducted by the “Overwhelmingly Large Simula¬ 
tions” (OWLS) project (Schaye et al. 2010) in a 100/i“^ Mpe box 
with 512^ dark-matter (and baryon) particles. 


Therefore, the mass resolution in our 512^ runs is compa- 
rable to that of K14, and our 7683 15363 resolutions 

are superior. Resolution can be important when applied 
to the large (100-200 kpc) Lya absorption systems ob¬ 
served in the low-redshift IGM. 

Star particles were formed in the simulations using the 
subgrid physics prescription of Gen & Ostriker (1992). 
Star formation occurs in cells where the baryon density is 
greater than 100 times the critical density, the divergence 
of the velocity is negative, and the cooling time is less 
than the dynamical time. As described by Smith et al. 
(2011), we adopted a star formation efficiency of 10% in 
these cells. Feedback from star formation returns gas, 
matter, and energy from star particles to the ISM and 
IGM. Although star particles are formed instantaneously, 
feedback occurs gradually with an exponential decay over 
time. Star particles produce 90% of their feedback within 
four dynamical times of their creation. Of the star par¬ 
ticles’ total mass, 25% is returned as gas, with 10% of 
that taking the form of metals and 10“® of the rest-mass 
energy returned as thermal feedback. The feedback is 
distributed evenly over 27 grid cells centered on the star 
particle. This scheme avoids the overcooling and subse¬ 
quent runaway star formation that occurs when feedback 
is returned to a single cell. We analyzed this “distributed 
feedback” method in our previous study (Smith et al. 
2011) and found that it produces better convergence and 
removes the over-cooling. These simulations were able 
to reproduce both the observed star-formation history 
and the number density of O VI absorbers per unit red- 
shift over the range 0 < z < 0.4. Many SPH studies of 
the IGM (Dave et al. 2010; Tepper-Garcia et al. 2012) 
suppress the over-cooling by suspending the cooling or 
turning off the hydrodynamics for a period of time, al¬ 
lowing the injected energy and metals to diffuse into the 
surrounding gas. 

The ionization states of H and He in the IGM were 
calculated with the non-equilibrium chemistry module in 
Enzo (Abel et al. 1997; Anninos et al. 1997). Metal cool¬ 
ing of the gas was computed using precompiled Cloud^ 
(Ferland 2013) tables that assume ionization equilibrium 
for the metals and are coupled to the chemistry solver 
(Smith et al. 2008). Photoionization and radiative heat¬ 
ing of the gas come from a spatially uniform metagalac- 
tic UVB. For each realization of the initial conditions, we 
ran three simulations that were identical except for the 
form of the UVB. Two simulations used analytic fits to 
the HMOl and HM05 backgrounds, initialized at z = 8.9 
and reaching full strength by z = 8. A third bacteround 
used the HM12 table implemented in the Grackl^chem- 
istry and cooling library (Bryan et al. 2014; Kim et al. 
2014) that begins at z = 15.13. 

In order to analyze the distribution of Lya absorbers in 
the simulations, we created synthetic quasar sight lines 
using the ytQ analysis package (Turk et al. 2011). For 
each simulation box, we created 200 synthetic sight lines, 
each constructed from simulation outputs spanning the 
redshift range 0.0 < z < 0.4. A randomly oriented ray 
was chosen through each output, representing the por¬ 
tion of the sight line beginning at the output redshift 

® http://nublado.org 

® https://grackle.readthedocs.org/ 

^ http://yt-project.org 
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and with sufficient pathlength to cover the Az to the 
next output hie. These rays were then combined to cre¬ 
ate a sight line spanning the entire redshift range. We 
used identical random seeds for simulations sharing the 
same initial conditions; differences in the absorber distri¬ 
butions therefore directly rehect effects of different ultra¬ 
violet backgrounds. Each sight line is composed of many 
line elements, from the portions of the sight line passing 
through a single grid cell. We identihed Lya absorbers 
as sets of contiguous line elements with neutral hydrogen 
density nni > cm“^. Egan et al. (2014) com¬ 

pared this method of identifying absorbers to a more so¬ 
phisticated method involving synthetic spectra. In their 
Figure 11, they showed that the resulting absorber col¬ 
umn density distributions are in good agreement. We 
also analyzed a large (1536^, Mpc) simulation run 

by Britton Smith and discussed in Shull et al. (2012a). 
The light rays from this simulation were generated by 
Devin Silvia, and we analyzed them in the same manner. 

2.4. Comparison of Simulations with Observations 

The flux decrement Da is an imperfect metric for es¬ 
timating the UVB, since measurements of the observed 
flux decrement depend on spectral resolution. In ad¬ 
dition, both measurements and simulations exhibit sig¬ 
nificant variance in Da- A more robust IGM diagnos¬ 
tic comes from the distribution of H I column densities. 
Figure 2 illustrates the observed STIS and COS distri¬ 
butions, /(IVhi, z), over the range 12.5 < logiVni < 14.5. 
The STIS and COS bins are offset by AlogAni = 0.1, 
using data from Tables 2 and 3. The COS survey is 
more extensive in both absorber numbers and column 
density, and its statistical accuracy is superior to that of 
the STIS survey. Over the overlapping range in log Ahi, 
the agreement is good. We over-plot values of /(A^hi, z) 
from our new IGM simulations (50/i“^ Mpc box and 
768^ grid cells) using three different ionizing backgrounds 
(HMOl, HM05, HM12). As analyzed further in Section 
2.5, the inferred baryon density requires ionization cor¬ 
rections for the neutral fraction, nui/u//. Appendix A 
suggests that the amplitude of the Lya absorber distribu¬ 
tion, /(Ahi, z), is proportional to for fixed 

baryon content (Df,) and temperature T = (10^ K)T 4 , 
and it therefore provides a constraint on the UVB. In 
our simulations, we find that /(Vhi,z) oc for ab¬ 

sorbers in the range 13 < log Vhi < 14. 

Figures 2-7 illustrate the differential absorber distribu¬ 
tion, dfN/d (log Nui)dz, with tests of convergence, cos¬ 
mic variance, box size, feedback, and redshift evolution. 
The simulated distributions are compared to observed 
distributions from HST surveys of Tilton et al. (2012) 
and Danforth et al. (2015). Figure 2 presents the pri¬ 
mary results of this paper: the H I column-density distri¬ 
butions for three radiation fields (HMOl, HM05, HM12) 
computed with 768^ grids and a 50h~^ Mpc box. The 
lower HM12 background produces significantly more ab¬ 
sorbers than earlier UVBs, ranging from a factor of two 
for absorbers between 12.5 < logAni < 12.7 to a fac¬ 
tor of seven for 14.3 < logAni < 14.5. The observed 
distributions fall between the results for the HMOl and 
HM12 backgrounds, approaching the HM12 simulations 
for log Ahi > 14.0. This suggests that there may be less 
tension between the HM12 background and observations 



log(N,j,) 


Fig. 2.— HST surveys of the column-density distribu¬ 
tion of low-redshift IGM absorbers are compared to new 
cosmological simulations with the Enzo grid code (768® 
cells, 50h~^ Mpc box) with three ionizing UV backgrounds 
(HMOl, HM05, HM12). We plot the bivariate distribution 
of absorbers, f{Nm,z) = d‘^jV/d(logNui)dz, in redshift and 
column density (Ahi in cm“®) with AlogAni = 0.2 bins. 
Data from HST/COS survey (D15: Danforth et al. 2015) 
are shown as red squares, and data from HST/STIS survey 
(T12: Tilton et al. 2012) as green crosses. Analytic theory 
(Appendix A) shows that, for fixed baryon content, Lya line 
frequency scales with photoionization rate and gas temper¬ 
ature as The HM12 simulations use a lower 

ionizing flux and produce a higher absorption-line density. 
For data with the best statistics (13 < log Ahi < 14), the ob¬ 
served distribution, /(Ahi,^) « (167)[Ahi/ 10®® cm“®]“°'®®, 
lies between simulations run with higher UVB (HMOl and 
HM05) and the lower HM12 background, except at the 
highest column densities. The best-fit distribution requires 
an ionizing background approximately 2-3 times higher than 
HM12, but with sufficient uncertainty that a crisis in photon 
production is unlikely. 


than implied by the K14 results. Figure 3 explores ef¬ 
fects of cosmic variance, comparing three models with 
50/i“® Mpc boxes and 512® grids. For each radiation 
field (HMOl, HM05, HM12) the simulations labeled with 
subscripts 1 and 2 are those with the Music initial con¬ 
ditions, and those labeled 3 use the Enzo-packaged initial 
conditions. For a fixed UVB the scatter among simula¬ 
tions is at the 15% level or less, indicating that variance 
is unlikely to play an important role. Figure 4 explores 
convergence of our simulations, comparing models with 
50/i-® Mpc boxes on grids of 512®, 768®, and 1536®. The 
differences between these models are relatively small, at 
the 10% level and within the expected differences arising 
from sample variance. Figure 5 investigates the potential 
effects of box size, which could suppress structure forma¬ 
tion down to low redshift if the box was too small. We 
ran two simulations (100/i“® Mpc, 512® and 50/i“® Mpc, 
256®), both using the same HMOl ionizing background. 
The grid cell sizes are identical (195h“® kpc), and the 
results are essentially the same for log Affii < 14. 

Figure 6 explores the potential influence of our feed¬ 
back methods. In our 512® simulations, feedback has 
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HM01_1 HM05_2 iTi: HM12_3 

HM05_1 HM12_2 ^ ^ D15 

IZZI HM12_1 iTr: HM01_3 i i T12 


: HM01_2 cr: HM05_3 

IqO I-^^^- 

12.5 13.0 13.5 14.0 14.5 

log(N^^) 



log(Ny^) 


Fig. 3.— Same format as Figure 2. To test cosmic variance, 
we show three (50/i“^ Mpc, 512^) simulations for each UVB, 
labeled HMOU, HM05_i, and HM12_i (subscripts i = 1,2,3 
correspond to runs with different initial conditions) and 
compared to HST survey data (D15 and T12). 



loglNa,) 


Fig. 4.— Same format as Figure 2. To test convergence, we 
compared three (50h“^ Mpc, 512^) simulations, labeled as 
HM05_i (subscripts i = 1, 2, 3) with two larger simulations: 
a new 768^ simulation (HM05_768) and the 1536^ simulation 
(HM05_1536) provided by Britton Smith and first presented 
in Shull et al. (2012). All used the same (HM05) ionizing 
radiation field and are compared to HST spectroscopic 
survey data from STIS (T12) and from COS (D15). 


little effect on the distribution of H I column densities. 
We find little difference if we inject mass, metals, and 
energy within a single cell, in adjacent cells, or elim¬ 
inating feedback entirely. Clumpiness of the IGM on 
small scales can be affected by the feedback method, even 
though the large-scale structural parameters, (T, pf,) and 
(Afe, Ahi), remain nearly the same. Figure 7 investigates 
the possible effects of redshift evolution of the Lya forest. 
We split the simulated results into two redshift intervals 
(0 < z < 0.2 and 0.2 < z < 0.4) in our 768^ simula¬ 
tions, all run with the HMOl radiation field. Differences 
between the simulated distributions are quite small and 


Fig. 5.— Same format as Figure 2. To test effects of 
box size, we show two simulations, one at (100/i“^ Mpc, 
512^) and another at (50/i“^ Mpc, 256®), both using the 
HMOl ionizing background. Results are compared to our 
HST survey data (D15 and T12). The grid cell sizes are 
identical (195h,“® kpc), and the results are essentially the 
same for log Nui < 14. The roll-off in simulated absorbers at 
logAni > 13.7 is an artifact of the poor resolution in these 
runs. 


. ^ i 

I— i 


i i ■ i 


■ i 


Distributed ^ ^ D15 
Locai i i T12 

None 


13.0 13.5 

loglNa,) 


Fig. 6.— Same format as Figure 2. To test effects of 
feedback injection methods, we show three (50h“® Mpc, 
512®) simulations following the feedback prescriptions of 
Smith et al. (2011). We explored distributed feedback (into 
the adjoining 27 cells), local feedback (into a single cell), and 
no feedback from star formation. All simulations are run 
with the same HMOl radiation field and compared to HST 
survey data (D15 and T12). 

consistent with COS observations (Danforth et al. 2015), 
in which the weak Lya absorbers (log IVhi < 14) exhibit 
little redshift evolution between z = 0 and z = 0.4. 

2.5. Comparisons to Previous Simulations 

To compare our grid-code simulations with previous 
SPH models, we examined two structural measures of the 
IGM. The first measure is the thermal-phase diagram of 
(T, pb), the temperature vs. baryon density in the low- 
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Fig. 7. — Same format as Figure 2. To test for redshift 
evolution of the low-redshift Lya forest, we split the results 
into two redshift intervals (0 < 2 < 0.2 and 0.2 < 2 < 0.4) for 
768® simulations, run with both HM05 and HM12 radiation 
fields. Evolution in the distribution of weak Lya absorbers 
in the HST survey (D15 and T12) are small, differing by less 
than 5% between 0 < 2 < 0.4. 

redshift Lya forest. For our grid models, these diagrams 
were shown as Figure 19 in Smith et al. (2011) and as 
Figure 1 in Shull et al. (2012b). For the SPH models, we 
examined Figure 8 of Dave et al. (2010) and Figure 7 of 
Tepper-Garcia et al. (2012). All of these diagrams are 
consistent with a relation T = (5000 K)A° ® where A;, is 
the baryon overdensity. A second measure is the correla¬ 
tion of baryon overdensity and column density (Af,, A^hi), 
typically expressed as a power-law, Af, = AqN^^, where 
Ag is the normalization at a fiducial H I column density 
A^hi = (10^"^ cm“^)Ai4. These correlations vary with 
redshift, but typically are evaluated at 2 = 0.25. Dave 
et al. (2010) used (384®,48/i“^ Mpc) simulations with 
the HMOl background to derive a fit for absorbers with 
T < 10'‘'5 K, 

Af, = (35.5 ± 0.03) . (6) 

Here, fr ~ 2/3 is a renormalization factor of opti¬ 
cal depths introduced to multiply the simulated op¬ 
tical depths to match the mean flux decrement Da- 
Thus, at 2 = 0.25 and with fr « 2/3, they found 
Af, = (38.9)Af°4^^®. From (512®, 100/i“^ Mpc) simula¬ 
tions with the HMOl background, Tepper-Garcia et al. 
(2012) found Af, = (48.3)with no renor¬ 
malization factor f-r- They noted a theoretical expecta¬ 
tion that Af, oc Afg/®® for absorbing gas in hydrostatic 
and photoionization equilibrium (Schaye 2001). In our 
(768®,50/i“® Mpc) simulations, also using the HMOl ra¬ 
diation field, we find Af, = (36.9)^°/® over the range 
0.2 < 2 < 0.3 for column densities 12.5 < log Ari < 14.5. 
The overdensity is calculated as a weighted average over 
the line elements that make up the absorber, using the 
column density as the weight field. We find similar nor¬ 
malizations, Ag = (36.9,38.9,48.3) for the three simula¬ 
tions (Shull et al. 2015, Dave et al. 2010, Tepper-Garcia 
et al. 2012) evaluated at 2 = 0.25 with the HMOl back¬ 
ground. 


Evidently, the structural measures of the low-2 Lya 
forest give similar results. The thermal-phase diagrams 
are essentially identical, while the differences in the 
(Af,, Ahi) correlation likely arise from different meth¬ 
ods of identfying and characterizing H I absorbers in 
column density and over-density. These comparisons 
are complicated by the fact that Dave et al. (2010) 
corrected their simulated optical depths by a factor 
(/r), and Tepper-Garcia et al. (2012) corrected for the 
temperature-density correlation of gas in hydrodynamic 
and photoionization equilibrium. The simulations also 
used different cosmological parameters. Tepper-Garcia 
et al. (2012) adopted density fractions (Dm = 0.238 and 
fib = 0.0418) from WMAP-3, which are lower by 15.6% 
and 9.1% than our WMAP-9 values (fl-m = 0.282 and 
fib = 0.046). Dave et al. (2010) adopted WMAP-7 pa¬ 
rameters (Dm = 0-28 and Df, = 0.046) similar to our 
values. Kollmeier et al. (2014) used Dm = 0.25 (10% 
lower than our value) and Df, = 0.044 (4.3% lower). The 
simulations also have different resolutions: our standard 
simulations use (768®, 50/i“® Mpc) while Tepper-Garcia 
use (512®, 100/1“^ Mpc). Given the scatter in these cor¬ 
relations, their redshift dependence, and the correction 
factors, the modest differences do not warrant concern. 

3. SUMMARY: A HIGHER UV BACKGROUND? 

Using our recent HST spectroscopic surveys of inter- 
galactic Lya absorbers, we have characterized their dis¬ 
tribution in H I column density and redshift. By com¬ 
paring the observed distribution to new Enzo simula¬ 
tions of the low-redshift IGM, we find an ionizing back¬ 
ground intermediate between the HMOl and HM12 val¬ 
ues. As shown in Appendix A, the inferred baryon den¬ 
sity and absorption-line frequency depend inversely on 
the UVB and temperature, scaling as Pjj®/^T-o.363^ 
fit the the HST data, we require approximately a fac¬ 
tor of 2-3 increase in the photoionization rate above 
HM12, somewhat less than the factor of five suggested 
by K14. However, no single UVB agrees with the full 
distribution. Figure 2 suggests that a higher UVB is 
needed to explain the line frequency of weak absorbers 
(12.7 < log A < 13.9), while the lower HM12 background 
is consistent with the distribution of stronger absorbers 
(logAni > 14). Given the uncertainties in source emis- 
sivities, cosmological radiative transfer, and galaxy LyC 
escape fractions, our H I results do not constitute a crisis 
in our understanding of the sources of the UVB. 

We have explored a number of potential explanations 
for the differences between the H I distributions pro¬ 
duced by various simulations and their sub-grid feedback 
schemes, gaseous sub-structure, and injection of mass 
and metals. In photoionized Lya absorbers, the neutral 
hydrogen density depends on Gonsequently, the 

UVB and photoionization rate needed to explain the 
amplitude of /(Ari, 2) may be sensitive to the clumping 
factor, Ch = However, a detailed code 

comparison is beyond the scope of what we can do 
at this time. In the comparisons discussed earlier, we 
demonstrated that the integrated column densities, 
Ari, are less sensitive to the procedures for identifying 
absorbers or the feedback schemes. After investigating 
convergence, cosmic variance, and feedback, we conclude 
that the differences must arise elsewhere. 
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The primary influence on the distribution of H I 
absorbers in the Lya forest is the ionizing radiation 
field. A larger UVB was also found in calculations that 
included contributions from both quasars and galaxies 
(Shull et al. 1999; Faucher-Giguere et al. 2009). The 
primary uncertainties in the UVB are the contribution 
from massive stars in galaxies (Topping & Shull 2015) 
and the LyC escape fraction (HM12; Benson et al. 2013), 
a highly directional quantity that is difficult to constrain 
statistically from direct observations. The parameteriza¬ 
tion adopted by HM12, /esc = (1.8 x 10“^)(1-|-z)^'^, was 
tuned to match observations at z > 2.5, but it is likely 
much too low at z < 2 where observational constraints 
are rare. Reliable values require direct measurements 
of LyC fluxes from a statistically significant sample 
(iVgai S> 20) of starburst galaxies to constrain the low 
escape fractions expected from a highly directional LyC 
escape geometry. For example, if theoretical models 
predict (/esc) ~ 0.05, with 5% of the LyC escaping 
from each side of a gaseous disk through perpendicular 
conical chimneys (Dove & Shull 1994), direct detections 
would be possible, on average, in only 5% of the obser¬ 
vations. Constraints on /esc will require large samples 
to deal with inclination bias. Such a situation is seen 
in direct-detection observations at z « 3 (Shapley et al. 
2006; Nestor et al. 2013; Mostardi et al. 2013) and at 
z « 0.2 (Heckman et al. 2011) which And either large 
(10-40%) fractions of transmitted LyC flux or none at 
all. The difficulties with low-z constraints on /esc were 
discussed by Shull et al. (2014) in their low-z census of 
IGM metal abundances and ionization ratios (C+^/C’*'^ 
and Si+3/Si+2). 

We therefore suggest that starburst galaxies are impor¬ 
tant contributors to the ionizing UVB at z < 2. Their 
contribution to the ionizing background would explain 
the observed distribution of Lya absorbers and could re¬ 
solve the discrepancy in IGM photoheating inferred from 
the opacity of the Lya forest at z < 2 (Puchwein et al. 
2015). It has also been suggested that TeV emission from 
blazars could add significant heat to the IGM through 
pair-production of high-energy electrons and positrons 
(Puchwein et al. 2012). We have not investigated their 
effects on our data, either on Da{z) or the distribution 
in H I column density. We summarize the main results 
of our study as follows: 

1. Compared to Enzo N-body hydrodynamical simu¬ 
lations of the low-z IGM with different values of 
the ionizing UV background, the observed distri¬ 
bution of Lya forest absorbers (z < 0.5) lies in¬ 
termediate between the HMOl and HM12 back¬ 
ground calculations. A fit to the observations re¬ 
quires a factor-of-two increase in the UVB above 
HM12 with a recommended hydrogen ionization 
rate rH(-2) = (4.6 x 10“^"^ s~^)(l + 

2. The one-sided ionizing photon flux $o ~ 
5700 cm“^ s“^ at z = 0 agrees with the ob¬ 
served IGM metal ionization ratios, C III/C IV and 
Si Ill/Si IV (Shull et al. 2014), and suggests a 25- 
30% contribution of Lya absorbers to the cosmic 
baryon inventory. 


3. The increased ionizing background probably re¬ 

quires an increase in the escape fraction of ionizing 
(LyC) radiation from starburst galaxies above the 
low values (/esc < 10“^) adopted by HM12. Ioniz¬ 
ing photons from galaxies with (/esc) ~ 0.05 gives 
results consistent with previous UVB modeling 
(Shull et al. 1999) that found similar contributions 
from AGN and starburst galaxies, with specific in¬ 
tensities (erg cm“^ s“^ Hz“^ sr“^) at 13.6 eV of 
Iagn = X 10-23 and /cai = Llto j x lO'^a. 

4. Because LyC propagation through the ISM is ex¬ 
pected to be highly directional, detections of es¬ 
caping photons depend on galaxy orientation and 
will require large surveys to obtain valid statistical 
inferences. 

Future observational and theoretical work could sig¬ 
nificantly improve our characterization of the ionizing 
UVB at z < 1. We need better accuracy of the column- 
density distribution of Lya absorbers at log Vhi < 13.0 
and log Vhi > 14.5. Larger surveys of IGM absorbers 
will allow us to measure the redshift evolution of the 
distribution, d^N'/d{logNHi)dz, to test the expected in¬ 
crease in photoionization rate, rH(2:) oc (1 -b z)^ "^. The 
most critical future experiment will be direct measure¬ 
ments of LyC escape fractions from a large sample of 
starburst galaxies at z < 0.4. 
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4. APPENDIX A: DEPENDENCE OF Lya ABSORBERS ON IONIZATION RATE 


The intergalactic Lya absorbers are thought to arise as fluctuations in dark-matter conflned clumps or filaments. 
Modeling their distribution and frequency throughout intergalactic space depends on the cosmological evolution of 
this gas, coupled to calculations of the “neutral fraction”, /hi = nm/nH, produced by the balance of photoionization 
and radiative recombination. These physical processes can be modeled by N-body hydrodynamic simulations, but one 
can also derive their general behavior through simple analytic arguments (Shull et al. 2012b). Here, we describe the 
key assumptions and parameters to illustrate how the baryon density and hydrogen ionization rate Th can be inferred 
from observations of the H I column density distribution and the frequency, dAf/dz, of absorbers per unit redshift. 

We assume a distribution of IGM absorbers with constant co-moving space density, (f){z) = /'o(l + z)^, and constant 
cross section . We adopt a cosmological relation between proper distance and redshift, dljdz = c{dt/dz) = 
[c/{l + z)H{z)], where the Hubble parameter at redshift z in a flat ACDM cosmology is H{z) = + 

In this model, the frequency of absorbers is dN/dz = [c/(l -I- z)H{z)] 4>{z). The internal density distributions of 

the absorbers are approximated as singular isothermal spheres, and the IGM baryon density and H I fractions follow 
from photoionization models that translate the observed H I to total hydrogen. To compare with observations, we 
integrate through the cloud along a chord at impact parameter p = (100 kpc)pioo of the AGN sight line through the 
absorber. We adopt a 100-kpc characteristic scale length of Lya absorbers, normalized at fiducial H I column density 
Nm = (lO^"* cm“^)A^i4. Because most Lya-forest absorbers are low density and optically thin in their Lyman lines, 

we use a case-A hydrogen recombination rate coefhcient, a^^ = (4.09 x 10“^^ cm^ , scaled to temperature 

T = (10^ K)T4. For nearly fully ionized H and He, with helium abundance y = riHe/nH ~ 1/12 by number, we adopt 
electron density Ue = {1 + 2j/)nH and mean baryon mass per hydrogen pi, = {1 + Ay)mH- 

The co-moving baryon mass density of Lya absorbers of column density A^hIj probed at impact parameter p, is the 
product of absorber space density. 


{dAf/dz)H{z){l +z) 
- 


(7) 


and absorber mass 


Mb(p) = 


2rH(z)jVHl(p) 

7r(l-F2i/)a^^ 


( 8 ) 


The co-moving mass density, pb = 4>oMb oc [p A^hi rHi]^^^(dA/’/dz), is then integrated over the distribution in H I 
column density. We normalize pb to the closure parameter, Q/ = pb/Per at redshift z = 0, relative to critical density 
Per = (SHq/SttG). Because our HST surveys extend to z « 0.5, we must include the cosmological evolution in H{z), 
(j){z), and rH(^). Before introducing these corrections, we outline the method at low redshift, for which the Euclidean 
calculation gives 


(HI) ^ 32(27r)^/^ 

^ 3 



(Gpb\ 

\cHo) 




1/2 


a^^\l + 2y) 


(7.1x10-3) 


dAf/dz 
100 


Pim A^i4 


Th 

rHM12 


1/2 

L —1 ^0.363 

AijO -^4 


(9) 

Here, we normalize the Lya frequency to a characteristic value dAf/dz = 100 and scale Th to the HM12 ionization rate, 
rHMi2 = (2.28 X lO”^"* s~^)(l + zY''^, which we fitted over the range 0 < z < 0.7. The weak temperature dependence 
comes from the recombination rate coefhcient. The coefhcient (7.1 x IQ-^) corresponds to 15% of the cosmological 
baryon fraction, Vtb = 0.0463 ± 0.00093 (Hinshaw et al. 2013; Ade et al. 2014). For fixed baryon density in the 

Lya forest, the observed frequency of Lya absorbers, dAf /dz oc F^ ' , scaling as the inverse square root of the ionizing 
background. 

Our HST surveys of Lya absorbers provide accurate values of based on an integration over the column-density 

distribution. After inserting the redshift corrections to (/(z), H{z), and rH(^), we derive a general expression for 0/ , 


n 


(HI) 


= (7.1 X 10-3) 


Eh 

1/2 

\H{z)] 

.rHM12. 


[ Ho \ 


1 1/2 


,-l„l/2 7^0.363 
Oo FlOO -‘-4 


a + z) 


0.2 


rN„ 


' 


d^Af 


d\og{Nm)dz 


TVif d(log/V hi) . (10) 


Owing to the rij^ dependence of neutral fraction in photoionization equilibrium, the factor rH(z) enters Equations (8) 
and (10) as the square root, with a (1 -I- z)^'^ dependence that nearly balances the (1 -I- z)^ from the cosmological ratio 
[(j){z)/{l-\-z)\ in Equation (7). Two recent HST surveys (Tilton et al. 2012; Danforth et al. 2015) probed Lya absorbers 
out to z « 0.47 and found that between 24% and 30% of the baryons reside in the low-z Lya forest. Equation (10) 
is consistent with those estimates, particularly with our recommended factor-of-two increase in ionizing background, 
rH(z) = (4.6 X IQ-^^ s~^)(l + zY''^ over the range 0 < z < 0.7. 



